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Abstract 

Background: Corals represent symbiotic meta-organisms that require harmonization among the coral animal, 
photosynthetic zooxanthellae and associated microbes to survive environmental stresses. We investigated 
integrated-responses among coral and zooxanthellae in the scleractinian coral Acropora formosa in response to an 
emerging marine pollutant, the munitions constituent, 1,3,5-trinitro-1,3,5 triazine (RDX; 5 day exposures to 0 
(control), 0.5, 0.9, 1.8, 3.7, and 7.2 mg/L, measured in seawater). 

Results: RDX accumulated readily in coral soft tissues with bioconcentration factors ranging from 1.1 to 1.5. 
Next-generation sequencing of a normalized meta-transcriptomic library developed for the eukaryotic components 
of the A. formosa coral holobiont was leveraged to conduct microarray-based global transcript expression analysis 
of integrated coral/zooxanthellae responses to the RDX exposure. Total differentially expressed transcripts (DET) 
increased with increasing RDX exposure concentrations as did the proportion of zooxanthellae DET relative to the 
coral animal. Transcriptional responses in the coral demonstrated higher sensitivity to RDX compared to zooxanthellae 
where increased expression of gene transcripts coding xenobiotic detoxification mechanisms (i.e. cytochrome 
P450 and UDP glucuronosyltransferase 2 family) were initiated at the lowest exposure concentration. Increased 
expression of these detoxification mechanisms was sustained at higher RDX concentrations as well as production 
of a physical barrier to exposure through a 40% increase in mucocyte density at the maximum RDX exposure. At 
and above the 1 .8 mg/L exposure concentration, DET coding for genes involved in central energy metabolism, 
including photosynthesis, glycolysis and electron-transport functions, were decreased in zooxanthellae although 
preliminary data indicated that zooxanthellae densities were not affected. In contrast, significantly increased 
transcript expression for genes involved in cellular energy production including glycolysis and electron-transport 
pathways was observed in the coral animal. 

Conclusions: Transcriptional network analysis for central energy metabolism demonstrated highly correlated 
responses to RDX among the coral animal and zooxanthellae indicative of potential compensatory responses to 
lost photosynthetic potential within the holobiont. These observations underscore the potential for complex 
integrated responses to RDX exposure among species comprising the coral holobiont and highlight the need to 
understand holobiont-species interactions to accurately assess pollutant impacts. 

Keywords: Coral holobiont, Marine pollution, Meta-transcriptomics, Acropora, Zooxanthellae, RDX, Symbiosis, 
Transcriptional network, Next generation sequencing 
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Background 

Coral reef ecosystems are experiencing dramatic world- 
wide declines due to the coupled impacts of local habitat 
degradation and global warming [1-3]. As these threats 
mount, the survival of corals increasingly depends on their 
capacity to rapidly establish physiological acclimation to 
environmental change [4-7]. While corals have shown re- 
markable resilience to extreme environmental changes 
through geological time [8], relatively little is known of the 
mechanisms that have facilitated acclimation. 

Environmental change in the form of pollutant release 
into marine environments also poses potential threats to 
coral survival. Historic and on-going military activities 
have resulted in the presence of discarded and unexploded 
ordinance (UXOs) containing the high explosive 1,3,5-tri- 
nitro-1,3,5 triazine (RDX) in critical coral reef habitats (i.e. 
Vieques, Puerto Rico [9] and Oahu, Hawaii [10]) posing 
unknown threats to already embattled ecosystems as the 
UXOs age and potentially leach RDX. Therefore, under- 
standing the response of coral to ordinance compound ex- 
posures such as RDX is critical for the conservation and 
preservation of these reef ecosystems. 

Corals represent a complex community called the "coral 
holobiont" composed of the coral animal, algal symbionts 
(zooxanthellae) within the coral tissues, and microbes 
inhabiting the coral mucus [11-14]. The photosynthetic 
activity of the zooxanthellae has been observed to contrib- 
ute as much as 95% of the carbon requirement for daily 
respiration and growth of the coral holobiont [15]. In re- 
turn, the corals provide a protective environment for the 
zooxanthellae in addition to essential nutrients derived 
from heterotrophic feeding [15,16] and fixation by associ- 
ated microbes [17]. 

As a critical step toward deciphering changes in coral 
holobiont interactions that reflect acclimation to envir- 
onmental stressors, we have completed the first broad- 
scale meta-transcriptomics study of the coral Acropora 
formosa. In marked contrast to previous sequencing ef- 
forts where coral genomic DNA was physically separated 
from that of zooxanthellae before analysis (i.e. [18-21]), 
we have characterized the protein-coding genes for the 
in situ holobiont including both coral and zooxanthellae 
(the metatranscriptome). This approach provided the 
ability to investigate global transcript expression within 
an intact symbiotic relationship between coral and zoo- 
xanthellae and thus enable the observation of interactive 
metabolism among species. A. formosa represents a useful 
coral model species that can be easily cultivated for use 
in controlled laboratory experimentation. Furthermore, 
A. formosa is ecologically and physiologically similar to 
other ecologically relevant Acropora species occurring 
in reef tracts around the world [22]. 

Here, we examined coral holobiont responses to RDX ex- 
posure using bioconcentration measures, coral-zooxanthellae 



meta-transcriptomics, transcriptional network analysis 
and preliminary studies of zooxanthellae and mucocyte 
density. In addition to providing meta-transcriptome 
characterization for the eukaryotic components of the 
A. formosa coral holobiont, this study demonstrates 
the integrative responses among the coral animal and 
zooxanthellae that comprise the coral holobiont to an 
emerging marine pollutant. 

Methods 

Ethics statement 

The work described in this paper represents laboratory 
studies conducted using aquarium-cultured Acropora 
formosa coral fragments that were purchased from the 
Oceans, Reefs and Aquaria Company (ORA, Harbor 
Branch Oceanographic Institution, Ft. Pierce, FL). The 
study did not involve vertebrate testing or experiments 
with threatened or endangered species. 

RDX exposures 

Fragments of the branched coral A. formosa ranging in 
size from 4 to 7 cm were allowed to acclimate to expos- 
ure chambers for 24 hours prior to experiment initiation 
(details on shipping and processing of corals is provided 
in the Additional file 1). Exposure chambers consisted of 
38 L glass aquaria filled with 20 L of reconstituted sea- 
water (Crystal Sea 8 Marinemix, Marine Enterprises Inter- 
national, Baltimore, MD) and equipped with a Coralife 
Super Skimmer and Bio Balls (Central Garden and Pet 
Co., Walnut Creek, CA) placed in the outlet filtration box. 
Given facility limitations and the requirement of this spe- 
cific exposure chamber setup to maintain coral health, the 
experimental design included single exposure chambers 
containing five independent biological replicates. While 
single exposure chambers do not allow exclusion of vari- 
ation in responses due to the exposure apparatus versus 
RDX exposure concentration, our principle concern was 
variation in biological responses at the individual level. 
Water quality, analytical chemistry and environmental 
conditions (methods and results described below) were 
closely monitored for consistency across all chambers to 
ensure variations in responses were due to the RDX ex- 
posure concentrations. We considered the individual coral 
fragments within each tank to represent true statistical 
replicates as has been described previously [23-25] only 
after careful empirical consideration through monitoring 
criteria contributing to statistical independence [26] and 
establishing dose-response relationships across the in- 
dependent exposure aquaria (see Additional file 1 for 
detailed discussion). 

Exposures included a control chamber (0 mg/L) and 
RDX-exposure chambers where RDX dissolved directly 
in seawater was added to exposure media targeting 0.5, 
1.0, 2.0, 4.0 and 8.0 mg/L concentrations (see Additional 
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file 1 for comments on expected exposure levels in the 
environment). Exposure media were maintained at 27°C 
by water re-circulating REMCOR heating/cooling units 
(REMCOR Products Company, Glendale Heights, IL). 
Salinity was maintained at 32 ppt and a 16:8 hr light-to- 
dark photoperiod applied with four high-intensity full- 
spectrum Phillips, Alto Collection F40T12/DX, 40 Watt, 
6500 Kelvin, 48 inch light bulbs (Phillips Corporation, 
Amsterdam, Netherlands) housed in high reflectance 
Lithonia Lighting 8 Model 1241DP chromed light fixtures 
(Acuity Brands Company, Conyers, GA) placed ap- 
proximately 7 cm above each exposure chamber. RDX 
concentrations and water quality parameters including 
temperature, salinity, pH, dissolved oxygen concentration, 
ammonia concentration, and nitrogen concentration were 
measured daily. At the termination of the RDX exposure 
bioassay, 3 of the 5 coral fragments from the control as 
well as each of the RDX exposure conditions were selected 
at random and immediately fixed in RNA Later™ (Ambion, 
Austin, TX) following manufacturers recommendations. 
For the remaining coral fragments, one was selected at 
random for RDX bioaccumulation analysis in coral soft 
tissues, while the final fragment was used for prelimin- 
ary investigations of RDX effects on coral soft tissue 
histology/histochemistry. 

Analytical chemistry of exposure water 

The concentration of RDX in exposure chambers was 
determined by daily sampling of exposure water. Two mL 
of exposure water was pipetted into 4 mL amber HPLC 
sampler vials and analyzed following USEPA method 8330 
[27] on a Waters HPLC using C18 and CN columns 
(detailed description provided in the Additional file 1). 
All RDX results are from the Supelco C18 column for 
which there were no interferences. The laboratory report- 
ing limit for all analytes was 0.5 umol/L (~0.1 mg/L) for 
water samples. 

Normalized cDNA library construction 

RNA extractions were conducted by gently scraping the 
coral soft tissue from the calcium carbonate skeleton and 
then following the Qiagen, RNeasy* Mini RNA extraction 
kit protocol (Qaigen Inc. Valencia, CA). The quantity and 
quality of the RNA extracted from each tissue sample was 
confirmed as described in Gust et al. [28]. The A. formosa 
RNA library was assembled by pooling all RNA samples 
from every replicate sampled in both control and all RDX 
exposure treatments. Reverse transcription of full length 
eukaryotic cDNAs was accomplished using Clontech 
SMART-PCR (Clonetech Laboratories Inc. Mountain 
View, CA). Due to the use of poly(d)T primers in this 
process, prokaryotic sequences which lack poly(A)-tails 
were excluded from the coral holobiont library. In order 
to generate sequence information for as many unique 



cDNAs as possible, the cDNA library was normalized 
using the Trimmer cDNA Normalization Kit (Evrogen 
JSC, Moscow, Russia) to insure sequencing of both high 
and low abundance transcripts. 

cDNA library sequencing and annotation 

The normalized cDNA library was sequenced using 
Roche 454 GS-FLX pyrosequencing [29] with several 
modifications that improved sequence reproducibility 
while reducing labor-intensive sample manipulations 
[30,31]. Flowgrams from the GS-FLX were assembled 
by Newbler Assembler (Roche Applied Science). Contigu- 
ous sequences (contigs) and singletons were searched 
against the NCBI Genbank protein database using blastx 
where E < 10 was considered a significant match to an 
expressed sequence tag (EST). Gene targets were anno- 
tated where possible with identity and functionality using 
Kyoto Encyclopedia of Genes and Genomes (KEGG, 
http://www.genome.jp/kegg/) and EuKaryotic Ortholo- 
gous Groups (KOG, http://genome.jgi.doe.gov/Tutorial/ 
tutorial/kog.html) databases. 

Phylogenetic characterization of the coral holobiont 
meta-transcriptome 

Phylogenetic designations connected to best blastx se- 
quence matches were generated for all contigs and single- 
tons having significant blastx scores (E<10~ b ). These 
phylogenetic identities were manually searched against 
published literature and assigned to one of three taxo- 
nomic groups (1) coral, (2) algal symbiont or (3) other. 
The group represented as "other" consisted of sequences 
where the phylogenetic origin was not representative of 
the two major eukaryotic groups contributing to the coral 
holobiont (i.e. fungi), or could not be definitively assigned 
to only one of these groups (i.e. various protists). All 
matches to prokaryotes were removed from the eukaryotic 
library. A "strict taxonomy" and an "inclusive taxonomy" 
were developed given our phylogenetic inquiry. The strict 
taxonomy only included associations beginning at the 
phylogenetic class anthozoa to represent coral and algal 
symbionts only included Symbiodinium species as well as 
algal species recognized to have direct associations with 
coral, all based on primary literature search. For the in- 
clusive taxonomy, all matches to animal species were 
included within the coral taxon and all plant and alga 
species were included within the algal symbiont taxon. 
Finally, given that the meta-transcriptomic library for 
the coral holobiont was normalized, we tested if the rela- 
tive expression of unique transcripts between the two 
eukaryotic contributors to the holobiont (coral and algal 
symbionts) was proportional using a binomial test (SAS 
v.9.2, SAS Inc. Raleigh, NC). 



Gust et al. BMC Genomics 2014, 15:591 
http://www.biomedcentral.com/1471-2164/15/591 



Page 4 of 18 



Development of custom microarray 

The 7,500 sequences identified with the lowest E values 
(which represent the best similarity to known genes) 
were uploaded to eArray (Agilent Technologies, Palo 
Alto, CA) in both forward and reverse sequence orienta- 
tion (due to lack of reading-frame information to assess 
the "sense" strand). eArray was utilized to incorporate 
these 15,000, 60mer oligonucleotide probes into a cus- 
tom microarray design (Agilent Technologies) for the 
eukaryotic members of the A. formosa holobiont. This 
design represents 7,500 sense-antisense probe-pairs pro- 
viding quality control of the microarray analysis for cross 
hybridization. Specifically, one probe in each probe pair 
represents a nonsense sequence and our expectation was 
that no target should have specific binding to it (see below 
for anti-sense strand exclusion methods). The phylogen- 
etic characterization of transcriptome sequence informa- 
tion (see previous section) was leveraged to provide 
taxonomic characterizations for all probes on the micro- 
array. The custom microarrays were developed for the A. 
formosa holobiont prior to the removal of sequences with 
blastx matches to prokaryotes, therefore 5.8% of micro- 
array probes represent non-quality assured sequence 
matches. Expression data for these probes was removed 
prior to final microarray analysis. 

Microarray hybridizations 

The effects of RDX exposure on transcript expression in 
the eukaryotic components of the holobiont (coral animal 
and algal symbiont) were assessed simultaneously by in- 
vestigating quality-assured RNA [32] collected from whole 
corals exposed in the exposure bioassay. The microarray 
hybridization experiment included 3 biological replicates 
for each of the following conditions: control, 0.5, 2, and 
8 mg/L, RDX (target concentrations). Microarray hybrid- 
izations were conducted utilizing the Agilent QuickAmp, 
One-Color Microarray Hybridization protocol (Agilent 
Technologies) following manufacturer's instructions. 

Microarray analysis 

An Axon GenePix* 4000B Microarray Scanner (Molecular 
Devices Inc., Sunnyvale, CA) was used to scan microarrays 
at 5 urn resolution. Data were extracted from microarray 
images using Agilent Feature Extraction software, version 
9.5.1 (Agilent Technologies). Microarray data were nor- 
malized to the 50th percentile within each array followed 
by median scaling among all exposures using GeneSpring 
Software version 7.3 (Agilent Technologies). GeneSpring 
was additionally used to conduct statistical analyses to 
identify differentially expressed transcripts (DET) among 
treatments using one-way ANOVA including Benjamini 
and Hochberg multiple testing corrections [33]. A 
post-hoc test including a parametric t-test (p = 0.05) 
and log 2 fold change cutoff of > 1.5 was used to discern 



statistically significant differences in transcript expres- 
sion for each RDX treatment relative to the control. 
Given the sense-antisense architecture of the microarray 
(see above), all differentially expressed microarray targets 
were examined to identify and eliminate microarray tar- 
gets for which both the sense and antisense probes were 
differentially expressed. Microarray data are publicly avail- 
able at the Gene Expression Omnibus (GEO, http://www. 
ncbi.nlm.nih.gov/geo/info/linking.html) under series ac- 
cession GSE27624. 

Delineating meta-transcriptomic expression in response 
to RDX exposure 

A Chi square test (Sigma Stat, v.3.1.1, Systat Software, 
Inc Chicago, IL) was used to test for differences in the 
proportion of transcripts differentially expressed among 
the phylogenetic groups (coral, algal symbiont and "other") 
across RDX-exposure concentrations. Chi-square tests 
were also used to determine if the number of differen- 
tially expressed transcripts changed in response to RDX 
concentration for each phylogenetic group. Similarly, 
we tested if the proportion of DETs changed among 
phylogenetic groups when investigating the number of 
DETs relative to the total number of probes represented 
for each group on the microarray. Finally, binomial tests 
(SAS v.9.2, SAS Inc. Raleigh, NC) were conducted to 
determine if the proportion of significant transcripts 
that had increased versus decreased expression deviated 
from the proportion expected by chance (i.e. 1:1 ratio). 

An investigation of transcript functions affected by RDX 
exposure among the phylogenetic groups was also con- 
ducted. Specifically, KOG, a eukaryotic representation of 
the Clusters of Orthologous Groups (COG) database [34], 
identities for all differentially expressed transcripts were 
leveraged to determine which metabolic pathways repre- 
sented the primary functional targets impacted by RDX 
exposure. 

Reverse engineering of transcriptional network 

We used an ordinary differential equations-based method 
described by Lai et al. [35] to construct an integrated tran- 
scriptional network. The network included all DET from 
both coral and zooxanthellae to assess correlations in ex- 
pression within and among species for the most highly- 
represented KOG functions. A visual representation of the 
resulting network was generated using Cytoscape [36]. 

Reverse-transcriptase, quantitative polymerase chain 
reaction (RT-qPCR) 

Microarray results were validated by RT-qPCR for 13 
unique transcript targets found to be of critical import- 
ance for the interpretation of RDX effects in the coral 
holobiont in addition to one regulatory "control" tran- 
script, "actin gene" for Acropora millepora (Additional 
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file 2: Table SI). Transcript expression levels were ex- 
amined using DNase (Qiagen, Valencia, CA) treated 
total RNA from each of the three biological replicates 
used in microarray hybridizations (see Additional file 1 
for methods). RT-qPCR data was analyzed with SDS 2.2 
software package (Applied Biosystems, Foster City, CA) 
using the AACT method to quantify results as recom- 
mended by the developer. Assumptions inherent in the 
AACT analysis were observed to be met sufficiently in- 
suring accuracy of downstream analysis. Actin is recog- 
nized to have relatively stable expression independent 
of external stimuli in cnidarians [37]. One-way ANOVA 
on threshold cycle (Ct) values including control, 0.5, 1.8 
and 7.2 mg/L exposures indicated that results are consist- 
ent with the assumption that actin expression remained 
consistent across treatments (p = 0.849), and was therefore 
selected as the internal normalizer for relative quantifica- 
tion (RQ). The characteristics of actin expression included 
low variability in relative quantification (RQ) and Ct value 
near the median of all RT-qPCR reactions. Fold change 
values (log 2 ) were calculated using RQ results where 
values represent transcript expression in RDX-treated 
coral relative to control coral. The 95% confidence 
interval (95% C.I.) was calculated around the mean 
relative expression for each RDX treatment. Confidence 
intervals that did not include unity were considered dif- 
ferentially expressed relative to controls as described in 
Rawat et al. [38]. 

Analytical chemistry of RDX-tissue residues 

One coral fragment per treatment was sampled for body 
residue determination. The soft coral tissues were gen- 
tly scraped from each A. formosa fragment with a sterile 
scalpel and transferred to QBiogene - Lysing Matrix A 
2.0 ml tubes (MP Biomedicals LLC) containing garnet 
matrix and one 0.6 cm ceramic sphere. HPLC-grade 
acetonitrile (0.25 ml, Sigma- Adrich) was added to each 
tube. The tissue samples were weighed and then homoge- 
nized in a FAST Prep-24 mini-bead-beater instrument 
(MP Biomedicals, Santa Ana, CA) for one minute and 
centrifuged at 7500 g for three minutes. After centrifu- 
gation, 0.1 mL of the extract was removed and mixed 
with 0.1 mL of CaCLi (5 g/L). Subsamples were added 
to 0.15 mL glass inserts and placed into 4 mL amber 
HPLC sampler vials. The samples were analyzed fol- 
lowing USEPA method 8330 [27] on a Waters HPLC 
using C18 and CN columns (detailed description pro- 
vided in the Additional file 1). All RDX results are from 
the Supelco C18 column for which there were no inter- 
ferences. Because RDX body residue was measured in 
only one replicate coral per treatment, statistical com- 
parisons were not conducted. The laboratory reporting 
limit for all analytes was 5 umol/kg (~1 mg/kg) for tis- 
sue samples. 



Preliminary study of coral tissue histology and 
histochemistry 

Changes in the coral tissue density of mucocytes and symbi- 
otic zooxanthellae in response to RDX exposure was investi- 
gated using ultra-high resolution analysis of coral holobiont 
histochemistry. Three technical replicate cuttings from a 
single coral fragment were collected from each exposure 
condition and fixed in Carnoys solution (60% ethanol, 30% 
chloroform and 10% glacial acetic acid). Methods for CaC0 3 
skeleton demineralization, tissue embedding, histological tis- 
sue sectioning, histochemical staining, image acquisition 
and analysis are described in Piggot et al. [39]. Briefly, histo- 
logical sections were individually imaged for fluorescence 
with mucocyte cells represented by AZ-acetylglucosamine 
content bound by wheat germ agglutinin-conjugated Alexa 
Fluor 647 and zooxanthellae cells represented by chloro- 
phyll auto-fluorescence induced with a metal halide bulb 
emitting through a HQ480/50 excitation filter. Mucocyte 
and zooxanthellae densities were quantified as described in 
Piggot et al. [39] with means and standard deviations calcu- 
lated for control and the 8 mg/L (target concentration) 
RDX exposure treatment. The standard deviation repre- 
sented the technical variation within a single replicate. 

Results 

Exposure conditions 

The concentration of RDX in the water remained rela- 
tively stable during the 5-d exposure (Additional file 2: 
Table S2a) where mean measured exposure concentrations 
were 0.50, 0.93, 1.77, 3.67 and 7.18 mg/L. Measured 
concentrations are used to represent RDX treatments 
throughout the remaining text. The greatest decrease in 
concentration between days 0 and 5 was 9%, observed 
in the 1 mg/L treatment. All water quality parameters 
were within adequate range according to guidance in 
ASTM [40], (Additional file 2: Table S2b). 

Coral meta-transcriptome sequencing, assembly and 
annotation 

The sequencing effort yielded 144.61 megabases in 702,750 
reads. After removal of adapter sequences, the assembled 
reads yielded 61,691 contiguous sequences (contigs) and 
127,925 singletons. Contigs and singletons were searched 
against the NCBI Genbank protein database using blastx 
identifying 12,141 significant gene matches at E < 10" 
after removal of microbial sequences (Additional file 2: 
Table S3). Gene targets were annotated where possible 
with identity, functionality, gene ontology, KEGG orthol- 
ogy and KEGG pathway information. 

Phylogenetic characterization of the coral 
meta-transcriptome 

Examination of the strict taxonomy generated for the 
coral meta-transcriptome indicated that more transcripts 
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could be definitively ascribed to coral as compared to 
algal symbionts (P < 0.001, Figure 1, Additional file 2: 
Table S3). The inclusive taxonomy provided a similar result 
where most transcripts were matched to coral (P < 0.001), 
however the inclusive taxonomy provided a broader repre- 
sentation of transcripts with lineages likely related to algae 
symbionts (Figure 1). The taxonomic distribution of tran- 
scripts represented as probes on the A. formosa holobiont 
microarray largely matched the distribution observed 
within the meta-transcriptome (Figure 1, Additional file 2: 
Table S4). The taxonomic distribution served as a refer- 
ence point for investigating the relative contribution of the 
coral animal and the algal symbionts regarding transcript 
expression discussed below. 

Although only 70 transcripts had best blastx matches 
to the Acropora genus, 6,943 best matches were made 
with Nematostella vectensis (Additional file 2: Table S3), a 
sea anemone with a fully sequenced genome [41] to which 
A. formosa is related at the phylogenetic class anthozoa. 
Matches to zooxanthellae included 10 sequences matching 
Amphidinium carterae and 56 matching species of the 
genus Symbiodinium which are the primary photosyn- 
thetic endosymbionts of Acroprid corals [42]. The major- 
ity of these matches were to the Symbiodinium clade 3. 
Expanding upon the strict taxonomy to an inclusive tax- 
onomy (Additional file 2: Table S3) allowed for vastly 



increased functional associations to the two predomin- 
ant eukaryotic classes of the coral holobiont community 
(corals and symbiotic algae) and provided a novel per- 
spective for understanding dynamics within the holo- 
biont in response to the RDX stressor exposure. 

RDX-tissue residues 

The concentration of RDX in coral tissues increased 
with increasing exposure concentration (Figure 2). The 
bioconcentration factor (BCF) values for RDX were rela- 
tively similar among treatments, with no apparent rela- 
tionship with the RDX concentration in water (Figure 2; 
mean BCF = 1.3). The mean BCF for A. formosa exposed 
to RDX was within the range of values observed for a 
variety of aquatic invertebrates and fish [43] . 

Microarray results overview 

As an initial quality control step, a total of 14, 4 and 36 
microarray targets within the 0.5, 1.8 and 7.2 mg/L RDX 
treatments, respectively, were identified to be differentially 
expressed in both the sense and antisense microarray 
probes when searched across all differentially expressed 
transcripts (DET). These targets represented non-specific 
probes and were therefore eliminated from the overall 
results set. The summation of the DET within each RDX 
treatment indicated a positive relationship between RDX- 
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Figure 1 Taxonomy of meta-transcriptome sequence data generated for the Acropora formosa coral animal and associated 
zooxanthellae. Taxonomies are provided for the meta-transcriptome as well as for the probes represented on the A formosa meta-transcriptome 
microarray. A strict and inclusive taxonomy are provided based on criteria explained in the Methods. 
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RDX BCF (L/kg) 
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RDX Target Exposure Cone. (mg/L) 

Figure 2 Measured RDX concentrations in exposure chambers 
and in coral tissue. Additionally, the bioconcentration factor (BCF) 
for RDX is provided. The y-axes for tissue residues and BCF are mg/kg 
and I7kg, respectively. 



exposure concentration and the number of DET (Figure 3) 
yielding a total of 120, 181 and 239 DET in the 0.5, 
1.8 and 7.2 mg/L RDX treatments, respectively (all 
differentially-expressed targets presented in Additional 
file 2: Table S5). The greatest degree of commonality in 
DET was found between the 1.8 and 7.2 mg/L RDX treat- 
ments (25.5%), followed by the 0.5 and 1.8 mg/L 



treatments (14.9%) whereas the least commonality (10.5%) 
was found among the most disparate RDX-exposure 
concentrations, 0.5 and 7.2 mg/L (Figure 3). Regression 
of log2 fold-change values for DET that were common 
among RDX treatments showed significant positive cor- 
relations (Figure 3). This indicates that within the DET 
there is a trend toward conservation of transcriptomic 
responses to RDX exposure across >1 order of magni- 
tude difference in RDX concentration. 

Differential impact of RDX among coral animal and 
zooxanthellae 

The proportion of total DET contributed among the coral 
animal and algal symbionts changed significandy across 
RDX-exposure treatments (p <0.001) where the relative 
contribution of algal symbionts increased with increasing 
RDX concentration (Figure 4). When normalized by the 
total number of probes represented for each taxonomic 
group on the microarray, the relative contribution of algae 
to DET was 58% and 57% in the 1.8 and 7.2 mg/L treat- 
ments, respectively (Figure 4). Multiple significant differ- 
ences in the proportion of transcripts having increased 
versus decreased expression were observed in response to 
RDX exposure concentration for coral and the algal sym- 
bionts (Figure 4). Coral had a greatly increased proportion 
of transcripts having increased expression at the 0.5 mg/L 
(p <0.001) and the 1.8 mg/L (p <0.001) RDX exposures, 
but not at 7.2 mg/L (p = 1.0). Finally, the proportion of 
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Figure 3 Effect of RDX on transcript expression in Acropora formosa. The bar chart and Venn diagram represent the total number of 
transcripts observed to have undergone significant differential expression compared to controls. The Venn diagram displays the number of 
differentially expressed transcripts that were common among the RDX treatments. Regression plots show log 2 fold-change relationships among 
differentially expressed transcripts found in common among RDX treatments. 
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Figure 4 Effect of RDX concentration on transcript expression in the coral and zooxanthellae. The pie chart labelled Total" provides the 
percent of the total suite of transcripts that had significant differential expression summed across coral, algal symbionts, and other ("other", 
representing unconfirmed taxonomic classifications). The pie chart labelled "Relative" provides the relative percentage of transcripts that had 
significant differential expression normalized to the total number of targets present for each taxonomic class on the microarray (See Additional 
file 2: Table S4 for target distribution on microarray). The bar graphs provide the proportion of significant transcripts that had increased vs 
decreased expression (totals inset in bars). Binomial tests were used to identify differences in the proportions of transcripts having increased : 
decreased expression with significant difference (p =0.05) marked with "*" and tests without sufficient data marked with "nt" representing no test. 



algal symbionts having decreased expression tended to in- 
crease with increasing RDX concentrations (Figure 4). 

Impacts of RDX on molecular pathways 

Investigation of KEGG functional terms for DET indi- 
cated the potential for impacts on a variety of molecular 
pathways in the RDX exposures (Table 1, and Additional 
file 2: Table S6). The predominant impact of RDX involved 
pathways underlying "metabolism". Specifically, carbohy- 
drate metabolism had the greatest number of total DET 
summing across all RDX exposure concentrations (Table 1). 
Differential expression of transcripts involved in amino 
acid metabolism was another prominent effect observed 



across RDX treatments. Additional prominent "metab- 
olism" pathways affected by RDX included: energy me- 
tabolism, lipid metabolism and metabolism of cofactors 
and vitamins. Aside from impacts on metabolism, impacts 
of RDX exposure on pathways involved in "environmental 
information processing" and "cellular processes" were add- 
itionally represented across all RDX treatments. 

Validation of microarray results using RT-qPCR 

Results of the RT-qPCR analysis corresponded directly 
with microarray results for 80.0%, 53.3% and 73.3% of 
targets for the 0.5, 1.8 and 7.2 mg/L RDX treatments, 
respectively (Figure 5). Further, significant positive 
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Table 1 Kyoto Encyclopedia of Genes and Genomes (KEGG) terms associated with transcripts that had significant 
differential expression in response to RDX exposure compared to controls in a 5d experiment 



0.5 mg/L RDX 
7. Metabolism 

1.1 Carbohydrate 
Metabolism 

Glycolysis/Gluconeogenesis 

Citrate cycle (TCA cycle) 

Pyruvate metabolism 

1.2 Energy Metabolism 

Nitrogen metabolism 

1.3 Lipid Metabolism 

Fatty acid metabolism 

Arachidonic acid 
metabolism 

1.4 Nucleotide Metabolism 

Purine metabolism 

7.5 Amino Acid Metabolism 

Glutamate metabolism 

1.6 Metabolism of Other 
Amino Acids 

Glutathione metabolism 



7.7 Glycan Biosynthesis and 
Metabolism 



Chondroitin sulfate 
biosynthesis 

Glycan structures - 
biosynthesis 1 



3. Environmental Information 
Processing 



Cor 
14 
5 



Alg Oth 1 .8 mg/L RDX (cont.) 



Pyruvate metabolism 



3. 7 Membrane Transport 7 0 

ABC transporters - General 
3.2 Signal Transduction 

MAPK signaling pathway 

ErbB signaling pathway 
Wnt signaling pathway 



Glyoxylate and dicarboxylate 
metabolism 

Butanoate metabolism 

7.2 Energy Metabolism 



0 0 Oxidative phosphorylation 

0 0 1.3 Lipid Metabolism 

0 0 Fatty acid metabolism 

0 0 Bile acid biosynthesis 

0 0 Glycerolipid metabolism 

0 0 7.5 Amino Acid Metabolism 

0 0 Glutamate metabolism 



0 0 Alanine and aspartate 
metabolism 

0 0 Glycine, serine and threonine 
metabolism 

0 0 Valine, leucine and isoleucine 
biosynthesis 

0 0 Arginine and proline 
metabolism 

0 0 1.6 Metabolism of Other 
Amino Acids 

0 0 Glutathione metabolism 



Cor Alg Oth 7.2 mg/L RDX (cont.) Cor Alg Oth 

0 0 1 1.4 Nucleotide Metabolism 1 7 0 

1 0 0 Purine metabolism 1 1 0 

0 0 1 7.5 Amino Acid Metabolism 4 0 3 

2 10 Alanine and aspartate 0 0 1 

metabolism 

2 1 0 Glycine, serine and threonine 2 0 1 
metabolism 

0 0 3 Valine, leucine and isoleucine 1 0 1 
biosynthesis 

0 0 1 Lysine degradation 



1 0 0 

10 0 2. Genetic Information 
Processing 

2.7 Transcription 

4 0 0 Aminoacyl-tRNA biosynthesis 

0 2.3 Folding, Sorting and 
Degradation 

0 0 Protein export 

0 0 

0 0 3. Environmental 

Information Processing 



1 0 0 3.2 Signal Transduction 

1 0 0 Two-component system - 
General 

Wnt signaling pathway 



0 0 1 7.6 Metabolism of Other 
Amino Acids 



0 0 1 One carbon pool by folate 



0 0 1 Porphyrin and chlorophyl 
metabolism 



0 0 1 Limonene and pinene 
degradation 



1 0 



1 0 0 Gamma- 

Hexachlorocyclohexane 
degradation 



1 0 0 

2 0 0 



0 0 1 Cyanoamino acid metabolism 1 0 0 
3 0 3 Glutathione metabolism 

2 0 0 7.9 Metabolism of Cof actors 
and Vitamins 



1 0 0 
7 7 0 

1 0 0 
0 



0 



1 0 0 

1 0 0 Phenylpropanoid biosynthesis 1 0 0 
Streptomycin biosynthesis 



1 0 0 
1 0 0 



Naphthalene and anthracene 1 0 0 
degradation 

1 Fluorene degradation 1 0 0 



0 0 1 

0 0 12. Genetic Information 302 
Processing 

7 0 0 2.1 Transcription 1 0 1 



1 0 0 Aminoacyl-tRNA biosynthesis 1 0 1 

2.2 Translation 2 0 1 

1 1 0 Ribosome 2 0 1 

7 7 0 

0 10 3. Environmental 9 1 0 

Information Processing 

1 0 0 3. 7 Membrane Transport 1 0 0 
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Table 1 Kyoto Encyclopedia of Genes and Genomes (KEGG) terms associated with transcripts that had significant 
differential expression in response to RDX exposure compared to controls in a 5d experiment (Continued) 
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KEGG terms are matched to source sequences: Coral (Cor), Algal Symbionts (Alg), or other (Oth) for non-specific phylogenetic associations. 



correlations were observed between transcript expres- 
sion among RT-qPCR and microarray results (Figure 5) 
indicating reasonable agreement among analytical tech- 
niques. Significant changes in expression were confirmed 
for a variety of DET representing a broad suite of mo- 
lecular functions from both coral and symbiotic algae 
(Figure 5). 



RDX effects on transcriptional networks 

The combined transcriptional network for coral and zoo- 
xanthellae generated using DET from the most promin- 
ently represented KOG terms, carbohydrate and energy 
metabolism, indicated several instances of highly corre- 
lated expression among the coral animal and the zoo- 
xanthellae in response to RDX exposure (Figure 6). 
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similar to UDPglucuronosyltransferase 2 family polypeptideB4 




0.57 
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Figure 5 Comparison of RT-qPCR and microarray results. Values represent tog2 fold change in transcript copy number relative to controls. 
Red and green highlighted cells represent statistically significant increases and decreases in copy number, respectively. Regression analyses 
represent correlations in log 2 fold change among microarray and RT-qPCR results. Linear regression significance tests, regression equations and R 2 
values provide evidence of correlations among the gene expression assays. "Species" represent putative source of transcripts. Primer sequences 
are provided in Additional file 2: Table SI. 



Given the intimate bioenergetic connection between 
coral and zooxanthellae where photosynthesis by zoo- 
xanthellae provides as much as 95% of energy resources 
to the coral holobiont [15], strong correlations in ex- 
pression within the combined network was expected. 
Examples of highly correlated expression among species 
included: pyruvate kinase (XP 001700637) in zooxan- 
thellae correlated with coral expression of transcripts 
involved in glycolysis (enolase 1 alpha, XP 001632906) and 
the electron transport chain (ATP synthase H + transporting 
mitochondrial Fl complex gamma, XP001638994 and 
cytochrome C oxidase subunit II, NP_6 12824). Within 
the complete transcriptional network, RDX exposure 
elicited strong interspecies network connections for 
genes involved in a broad compliment of energy metab- 
olism including: photosynthesis, glycolysis, pyruvate 
metabolism, the citric acid cycle, electron transport and 
oxidative phosphorylation (Figure 6, Additional file 2: 
Table S6). Although these results do not imply direct 
metabolic connections among these processes, the over- 
all trend of the results indicates a dichotomous response 
among species where expression of gene transcripts related 
to nearly all aspects of energy metabolism are decreased 



in zooxanthellae which is correlated with a broad scale in- 
crease in expression of energy metabolism-related gene 
transcripts in the coral. 

Preliminary study of coral tissue histology and 
histochemistry 

A. formosa exposed to RDX displayed no gross-level visual 
signs of pigmentation loss or excessive secretion of mucus. 
Comparison of mucocyte densities between the control 
(4720 ± 320 cells/cm 2 , mean ± SE) and coral exposed to 
7.2 mg/L RDX (7800 ± 1300 cells/cm 2 , mean ± SE) dem- 
onstrated the potential for increased mucocyte density at 
the highest RDX exposure thereby suggesting a possible 
mucosal protective response by the coral in response to 
RDX exposure (Figure 7). In contrast, histological examin- 
ation indicated no observable difference in coral tissue 
zooxanthellae densities between the control and samples 
exposed to 7.2 mg/L RDX as the means were similar with 
relatively little variance in either condition (Figure 7). 

Discussion 

The normalized cDNA library for A. formosa included 
transcripts (protein-coding RNAs) for two prominent 
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Figure 6 Transcriptional network inference of coral holobiont interactions. The network analysis demonstrates correlations of expression 
among molecular targets and species comprising the coral holobiont in response to RDX exposure. The network includes all Kyoto Encyclopedia 
of Genes and Genomes terms involved in carbohydrate and energy metabolism. Edge thickness represents correlation strength among nodes. 
Arrows represent "activation" from one node to the next while "T" ends represent "inhibition". At least one RDX exposure concentration elicited 
significant differential transcript expression for each gene, however we have provided expression information for all treatment levels to show 
trends. 



eukaryotic components of the coral holobiont (the coral 
animal and symbiotic zooxanthellae), thus representing a 
meta-transcriptome. Previous large-scale sequencing ef- 
forts for coral {Acropom millepora, [18,20] and Acropora 
digitifera [21]) and the close phylogenetic relative sea 
anemone, Nematostella vectensis, [41] were each sequenced 
without their associated algal symbionts. The inclusion of 
algal symbionts in our sequencing effort provided the 
ability to investigate the relative contribution of the coral 
animal and zooxanthellae to the meta-transcriptome as 
well as integrated transcript expression. 



The coral holobiont meta-transcriptome 

Both the strict and inclusive taxonomies generated for 
the coral holobiont meta-transcriptome indicated that 
the overwhelming majority of transcripts originated 
from coral relative to algal symbionts (Figure 1). We 
sought to determine if this observation was a result of 
differences in the number of protein coding genes in- 
herent in the respective genomes of each species. The 
closest phylogenetic relative of A. formosa having a 
fully sequenced genome, Acropora digitifera, was ob- 
served to have 23,700 protein coding transcripts [21]. 
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Figure 7 The effect of RDX exposure on zooxanthellae and 

mucocyte density in Acropora formosa. The bar chart shows the 

quantification of fluorescence from polyp histology sections 

showing auto-fluorescent zooxanthellae (green) and WGA-induced 

fluorescent mucocytes (red), 
k ) 



The recently sequenced transcriptome for the primary 
endosymbiotic dinoflagellate found in corals, Symbiodi- 
nium sp. reported 40,000-90,000 protein coding genes [44]. 
These observations compared to our meta-transcriptome 
results (Figure 1) indicate that contribution to the meta- 
transcriptome is not proportional to individual species tran- 
scriptome size, at least when comparing the eukaryotic con- 
tributors to the holobiont. This observation presents 
questions about individual-species transcriptomic expres- 
sion as it relates to holobiont function. For example, given 
the limited percentage of genes attributed to zooxanthellae 
in the total sequencing pool, is the diversity of transcript 
expression reduced in zooxanthellae engaged in symbiotic 
relationships compared to free-living zooxanthellae? Given 
that the zooxanthellae are engaged in a symbiotic associ- 
ation within the coral animal, division of physiological 
labor among species may have contributed to a decreased 
proportional expression of genes within zooxanthellae. For 
example, if essential nutrients required in zooxanthellae 
are directly supplied by coral, the need to express the 
molecular machinery innate in zooxanthellae to capture 
or create these metabolites is negated, thus the need to 
express the genes in such pathways is reduced. 

RDX bioaccumulates readily into coral tissue 

RDX bioaccumulated into coral tissue in a concentration- 
dependent manner and exhibited bioconcentration factors 
(BCF) ranging from 1.09-1.50 L/kg among RDX-exposure 
concentrations (Figure 2). A BCF near 1.0 L/kg indicates 
that RDX accumulation in whole-tissue samples mirrors 
ambient water concentrations, which is consistent with 
previous observations for RDX [43]. These results indicate 
that RDX was accumulated in coral tissue and was there- 
fore present to elicit potential toxicological affects within 
both the coral animal, as well as symbiotic zooxanthellae. 

Differential sensitivity among members of the holobiont 
at low RDX level 

Expression of transcripts in the zooxanthellae was relatively 
unaffected at the lowest RDX-exposure concentration 
(0.5 mg/L), whereas a relatively large number of transcripts 
were differentially expressed in the coral animal (Figure 4). 
At this concentration, RDX predominantly elicited differ- 
ential expression of various components of metabolism in 
the coral animal including carbohydrate, energy, lipid, 
amino acid and glycan metabolism (Table 1) as well as in- 
creased expression of potential detoxification mechanisms 
(i.e. cytochrome P450, XP_001624662, and UDP glucuro- 
nosyltransferase 2 family, XP 001638304, Figure 5). 

Carbohydrate metabolism was the most represented 
second order KEGG term in the 0.5 mg/L treatment. 
All transcripts involved in carbohydrate metabolism 
had increased expression (Additional file 2: Table S6) 
including phosphoenolpyruvate carboxykinase 1 (PCK1, 
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XP_001179448) which was RT-qPCR-validated (Figure 5). 
PCK1 has been observed to be the regulatory control 
point in the equilibrium reaction for glycolysis versus glu- 
coneogenesis [45] in a mammalian model, and if this func- 
tion is conserved in coral, the observed increase in PCK1 
expression indicates inertia toward glycolysis. Assuming 
orthologous gene functions in corals to those observed in 
the literature, the transcript expression results suggest an 
equilibrium shift toward increased energy consumption in 
coral in response to low concentrations of RDX. 

Cytochrome P450 metabolism is a key facilitator for 
chemical detoxification via phase I xenobiotic biotrans- 
formation [46] and UDP glucuronosyltransferase 2 rep- 
resents a xenobiotic phase II metabolism enzyme that 
conjugates sugar moieties to drugs and endogenous 
compounds increasing solubility and potential for elim- 
ination via excretion [47]. The qPCR validated increase 
in cytochrome P450 and UDP glucuronosyltransferase 2 
family transcriptexpression at this lowest RDX concentra- 
tion (Figure 5) indicates a likely response to metabolize 
bioaccumulated RDX. Although the bioaccumulation 
results suggest these mechanisms were not effective in 
decreasing RDX body burdens, such responses to chemi- 
cals where these detoxification mechanisms metabolize 
and mitigate toxicity would benefit both coral and the 
symbiotic zooxanthellae. Regarding the discussion of 
energy expenditures above, any induced detoxification 
mechanism requires energy which given the current ob- 
servations may account for some component of the in- 
creased inertia toward energy consumption at the lowest 
RDX exposure concentration. 

Integrated coral holobiont responses at elevated RDX levels 

Similar to the 0.5 mg/L exposure, effects on transcripts 
coding for genes involved in metabolism were the most 
pervasive responses at the 1.8 and 7.2 mg/L exposure levels 
(Table 1). However, in contrast to the coral-dominated re- 
sponse observed at the 0.5 mg/L treatment, the relative 
proportion of differentially expressed transcripts shifted to 
a zooxanthellae-dominated response with increasing RDX 
concentration (Figure 4). The most strongly represented 
metabolic function was carbohydrate metabolism in each 
the 1.8 and 7.2 mg/L RDX exposures, with additional 
marked representation of transcripts involved in energy 
and amino acid metabolism (Table 1). Cofactor/vitamin 
metabolism was additionally highly represented in the 
7.2 mg/L treatment. Each of these metabolic functions is 
discussed regarding the integrative response by the coral 
animals and zooxanthellae in the following. 

Carbohydrate metabolism 

Carbonic anhydrase activity has been observed to be re- 
quired by both species in invertebrate-algal symbiosis to 
support carbon fixation from C0 2 via photosynthesis 



[16] as well as playing a key role in coral skeleton depos- 
ition [48]. Transcripts coding for carbonic anhydrase had 
decreased expression at 7.2 mg/L in coral (AAM94169) 
and in zooxanthellae (AAW79301, Additional file 2: 
Table S5) indicating both reduced potential for coral 
skeleton growth and reduced potential to generate carbon 
substrates for metabolism via carbohydrate metabolic 
pathways. Regarding impacts on carbohydrate metabolism, 
decreased metabolic potential was also demonstrated in 
zooxanthellae where RT-qPCR validated reductions in 
pyruvate kinase (XP_001700637), the enzyme that cata- 
lyzes the reaction where phospoenolpyruvate is metab- 
olized yielding pyruvate and ATP [49] was observed 
(Figure 5, Additional file 2: Table S6). These results are 
suggestive of reduced potential for carbon fixation and 
cellular energy production via carbohydrate metabol- 
ism in the zooxanthellae as well as reduced potential 
for coral skeleton deposition in the RDX exposure. 

Changes in transcript expression related to carbohydrate 
metabolic pathways in the coral animal also reflected the 
potential for diminished substrate availability at the 1.8 
and 7.2 mg/L RDX exposure concentrations. First, expres- 
sion of transcripts coding for genes that metabolize large 
6-carbon carbohydrate substrates was reduced in coral. 
For example, RT-qPCR validated reduced expression of 
UDP-glucose 6-dehydrogenase (XP 971267) was observed 
with increasing RDX-exposure concentrations (Figure 5, 
Additional file 2: Table S6). UDP-glucose 6-dehydrogenase 
catalyzes reactions oxidizing the 6 carbon glucose sub- 
strate thereby reducing NAD + to NADH [50] which is a 
key energy-carrying molecule involved in cellular energy 
(ATP) production. Additionally, decreased expression of a 
gene transcript involved in metabolism of a 6 carbon 
intermediate substrate in the pentose phosphate pathway, 
6-phosphogluconate dehydrogenase (XP 001627413, [51]) 
was observed in the coral animal. In contrast to expression 
changes related to large 6 carbon carbohydrate substrates, 
expression of transcripts coding for genes that metabolize 
small 2- and 3 -carbon substrates had increased expression 
in corals. For example, at the 1.8 mg/L exposure concen- 
tration, citrate synthase (XP_001641037), a highly con- 
served enzyme known to regulate the first reaction in the 
citric acid cycle converting the 2 carbon substrate acetyl 
CoA to citrate [52] had increased expression in coral 
(Additional file 2: Table S6). Additionally, the coral animal 
had increased transcript expression of enolase 1 alpha 
(XP 001632906), the enzyme that catalyzes the reaction 
converting the 3-carbon substrate phosphoglycerate to 
phospoenolpyruvate at the 7.2 mg/L exposure concentra- 
tion. These results suggest that in the face of reduced car- 
bon fixation and energy production resulting from RDX 
exposure in zooxanthellae, larger carbohydrate substrates 
(i.e. glucose and pentose) may be in short supply causing 
corals to increase expression of enzymes that can capture 



Gust et al. BMC Genomics 2014, 15:591 
http://www.biomedcentral.com/1471-2164/15/591 



Page 15 of 18 



energy from intermediate 2 and 3 carbon metabolites to 
maintain energy budgets. 

Energy metabolism 

The observation that RDX exposure may cause diminished 
carbon fixation as indicated through reduced transcript 
expression of carbonic anhydrase in both coral and zoo- 
xanthellae (described above), it is intuitive that additional 
impacts on energy metabolism in the holobiont were ob- 
served. Reduced expression of transcripts coding photo- 
system II protein (CAI58622, Additional file 2: Table S6) 
as well as RT-qPCR validated reductions in pyruvate kin- 
ase (XP_001700637), chloroplast light harvesting complex 
protein (AAW79366) and major facilitator superfamily 
(AAU87579) in algal symbionts (Figure 5) are additional 
indicators of reduced potential for photosynthesis. 
Highly-stressed corals (particularly heat stressed) may 
eject zooxanthellae resulting in coral bleaching events 
[4] however, RDX exposure did not affect the density of 
algal symbionts (zooxanthallae) in coral tissue even at 
the highest RDX exposure concentration (Figure 7). 
Therefore, the decrease in transcriptional expression for 
zooxanthellae-associated genes, including those involved 
in photosynthesis, was unlikely an artifact of a bleaching 
event. Photosynthesis by zooxanthellae supports as much 
as 95% of energy metabolism in the coral holobiont [15]. 
Increased expression of gene transcripts involved in elec- 
tron transport chain-based ATP production were observed 
in the coral animal at both the 1.8 and the 7.2 mg/L RDX 
treatments likely representing a compensatory response to 
decreased cellular energy availability stemming from ap- 
parent decreased photosynthetic potential in zooxanthel- 
lae (Figure 6 and Additional file 2: Table S6). Given these 
expression results, we hypothesize that RDX may impair 
photosynthetic activity in zooxanthellae thereby affecting 
energy availability for the overall coral holobiont. In future 
studies, measuring photosynthetic rates in the holobiont 
under RDX exposures could explicitly test this hypothesis. 

Amino acid metabolism 

Changes in expression of transcripts involved in amino 
acid metabolism occurred in the coral animal, but not in 
zooxanthellae. RT-qPCR validated an inverse relationship 
between 3-methyl-2-oxobutanoatedehydrogenase lipoa- 
mide (XP 647496) transcript expression with increasing 
RDX-exposure concentration indicating reduced process- 
ing of amino acid substrates [53] and thereby reduced 
ATP consumption to fuel this process. These results are 
consistent with a metabolic shift to conserve energy 
resources. 

Metabolism of cofactors and vitamins 

Similar to observations at the lowest RDX exposure con- 
centration, increased expression of cytochrome P450 and 



UDP glucuronosyltransferase 2 family transcripts were ob- 
served in the RDX-exposed coral animal and validated 
using RT-qPCR (Figure 5 and Additional file 2: Table S6). 
In contrast, no transcriptional effects on cytochrome 
P450, UDP glucuronosyltransferase 2 family, or other 
phase I or II detoxification mechanisms were observed 
in zooxanthellae suggesting coral may be the principle 
component of the holobiont that actively initiated mecha- 
nisms to metabolize and potentially detoxify RDX. 

Mucus production in response to RDX 

The preliminary histochemistry assessment suggested that 
RDX exposure may have the potential to increase coral- 
derived mucocyte density at the highest RDX-exposure 
concentration (Figure 7). Corals have been observed to re- 
spond to a wide variety of environmental stressors such as 
crude oil, copper sulphate, mercury, increased sea surface 
temperature and decreased salinity by excreting copious 
amounts of mucus (see review in Brown and Bythell [13]). 
When stimulated, coral tissue mucocytes discharge hydro- 
scopic mucin into the seawater directly overlying the 
coral, creating an instantaneous large-volume expansion 
of mucus over a large surface area of coral tissue. This 
coral surface mucus layer has been observed to serve as 
a medium for nutrient assimilation, and is hypothesized 
to physically protect coral from chemical exposure by 
increasing the thickness of uncontaminated mucus in 
direct contact with the coral [13]. These preliminary re- 
sults in conjunction with the observed increases in gene 
transcript expression for xenobiotic detoxification mecha- 
nisms are suggestive of both a physical and metabolic re- 
sponse by the coral animal to protect the holobiont in a 
chemical stressor exposure. 

A hypothetical scenario for RDX-induced energy loss in 
the holobiont 

The critical role of zooxanthellae photosynthesis in gener- 
ating energy within the coral holobiont is well described 
[15,54]. Preliminary results indicate that holobiont bio- 
accumulation of RDX (Figure 2) did not affect zooxanthel- 
lae density in the A. formosa holobiont (Figure 7). RDX 
did however decrease transcript expression for genes in- 
volved in carbon fixation from CO2 (carbonic anhy- 
drase, AAM94169 and AAW79301, Additional file 2: 
Table S5), photosynthesis [photosystem II protein 
(CAI58622, Additional file 2: Table S6), chloroplast 
light harvesting complex protein (AAW79366, Figure 5) 
and major facilitator superfamily (AAU87579, Figure 5)], 
and carbohydrate metabolism (pyruvate kinase, XP_ 
001700637, Figure 5). Consistent with a decreased 
potential for photosynthetic energy production by zoo- 
xanthellae, coral animals decreased transcriptional ex- 
pression of genes involved large 6 carbon carbohydrate 
metabolism (UDP-glucose 6-dehydrogenase, XP 971267, 
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Figure 5 and 6-phosphogluconate dehydrogenase XP_ 
001627413, Additional file 2: Table S6) in favor of in- 
creased expression for metabolism of smaller 2 and 3 
carbon substrates in glycolysis (phosphoenolpyruvate 
carboxykinase 1, XP 001179448, Figure 5 and enolase 
1 alpha, XP_001632906, Additional file 2: Table S6) 
and the citric acid cycle (citrate synthase, XP 001641037, 
Additional file 2: Table S6). The transcriptional network 
generated for combined coral-zooxanthellae carbohydrate 
and energy metabolism showed strong correlations be- 
tween the decreased expression of photosynthetic genes in 
zooxanthellae and increased transcript expression for 
genes involved in glycolysis, the citric acid cycle and the 
electron transport chain (Figure 6). Although it is unlikely 
that these correlative network connections represent ca- 
nonical feedback in metabolism among coral and zooxan- 
thellae, they do indicate that the coral animal has highly 
sensitive and integrated transcriptional responses to the 
energetic potential of the zooxanthellae. 

Conclusions 

Our results indicate differential sensitivities and re- 
sponses to the emerging marine pollutant RDX among 
species comprising the A. formosa coral holobiont. The 
coral animal had the most sensitive transcriptional re- 
sponse to RDX bioaccumulation (Figure 4) including in- 
creased expression of genes involved in xenobiotic 
metabolism such as phase I and II xenobiotic detoxifica- 
tion pathways, cytochrome P450s and UDP glucurono- 
syltransferase 2 family enzyme, respectively (Figure 5) at 
the lowest RDX exposure level. These results in addition 
to preliminary data suggesting increased mucus produc- 
tion (Figure 7) as a physical barrier to RDX exposure indi- 
cate that the coral animal may provide a sentinel response 
to protect the coral holobiont against chemical stressors. 
As RDX exposure concentrations/bioaccumulation in- 
creased, the proportion of differentially expressed tran- 
scripts from zooxanthellae dominated the holobiont 
transcriptional response (Figure 4). Transcriptional 
network (Figure 6) and metabolic pathway analyses 
(Additional file 2: Table S6) suggest that the coral holo- 
biont may have reduced potential for carbon fixation and 
energy production via photosynthesis in zooxanthellae. 
Observed compensatory increases in expression of glyco- 
lytic, citric acid cycle and electron transport pathways 
were observed likely to sustain energy (ATP) needs in the 
coral animal. These observations underscore the potential 
for complex integrated responses that occur among spe- 
cies comprising the coral holobiont and highlight the need 
to understand holobiont-species interactions to accurately 
ssess stressor impacts. 
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